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Abstract: The conventional zenith tropospheric delay (ZTD) model (known as the Saastamoinen model) does not consider 
seasonal variations affecting the delay, giving it low accuracy and stability. This may be improved with adjustments to 
account for annual and semi-annual variations. This method uses ZTD data provided by the Global Geodetic Observing 
System to analyze seasonal variations in the bias of the Saastamoinen model in Asia, and then constructs a model with 
seasonal variation corrections, denoted as SSA. To overcome the dependence of the model on in-situ meteorological 
parameters, the SSA+GPT3 model is formed by combining the SSA and GPT3 (global pressure-temperature) models. The 
results show that the introduction of annual and semi-annual variations can substantially improve the Saastamoinen model, 
yielding small and time-stable variations in bias and root mean square (RMS). In summer and autumn, the bias and RMS 
are noticeably smaller than those from the Saastamoinen model. In addition, the SSA model performs better in low-latitude 
and low-altitude areas, and bias and RMS decease with the increase of latitude or altitude. The prediction accuracy of the 
SSA model is also evaluated for external consistency. The results show that the accuracy of the SSA model (bias: —0.38 cm, 
RMS: 4.43 cm) is better than that of the Saastamoinen model (bias: 1.45 cm, RMS: 5.16 cm). The proposed method has 


strong applicability and can therefore be used for predictive ZTD correction across Asia. 
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1. INTRODUCTION 


Tropospheric delay is the main error source for Global 
Navigation Satellite System (GNSS) positioning, navigation 
and timing. The effect of the troposphere on the GNSS signals 
appears as an extra delay in the measurement of the signal 
traveling from the satellite to receiver. The error caused by 
the troposphere is about 2 m in the zenith direction and 
20 m for lower elevations!!], and this must be corrected in 
GNSS applications. Unlike ionospheric delay!!, tropospheric 
delay is not dependent on signal frequency, and consequently 
cannot be eliminated by combining signals of different fre- 
quencies. Therefore, an empirical model of zenith tropo- 
spheric delay (ZTD) is usually used to correct it. Accord- 
ingly, it is of great practical significance to construct a high- 
precision ZTD empirical model. At present, the commonly 
used ZTD empirical models can be divided into the following 
two categories: 

(1) Meteorological parameter models. These models 
use meteorological parameters such as surface pressure, tem- 
perature and water vapor pressure to calculate ZTD. Among 
these models, the accuracy of the Hopfield model decreases 


with increasing altitude, while the Saastamoinen model is 
affected little by altitude, so the latter is more widely applica- 
ble. The ZTD correction accuracy of these models can 
approach centimeter precision with measured meteorological 
parameters, but the correction accuracy can decrease 2 to 3 
times if the standard meteorological parameter model is 
used, which can limit the applicability of these models in 
high-precision GNSS positioning, navigation and timing. 
To remove dependence on measured parameters, Liu et al., 
Yang et al., Yao et al., and Du et al. use the meteorological 
parameters provided by the global barometric temperature 
GPT/GPT2/GPT2w model as the input of the Saastamoinen 
modell5-¢l], which effectively improves the applicability of 
these two models in high-precision GNSS applications. 

(2) Non-meteorological parameter models, such as the 
European Geostationary Navigation Overlay Service 
(EGNOS)!’: 81, which are the ZTD correction model adopted 
by the European and American wide-area augmentation sys- 
tems. They only need station location information, without 
surface meteorological parameters, to calculate ZTD, and 
the global average ZTD calculation accuracy is comparable 
to that of Saastamoinen and Hopfield models based on mea- 


sured parameters. However, the EGNOS and University of 
New Brunswick (UNB) models divide the earth into 15° lati- 
tude intervals, so the spatial resolution is poor and the effect 
of longitude is ignored, resulting in an inability to capture 
local variation in ZTD. Furthermore, the ZTD correction cal- 
culated by these models is ineffective in areas outside 
Europe and North America. To meet the demand for ZTD cor- 
relationin high-precision GNSS applications, new non-meteo- 
rological parametric ZTD models have been successively 
established. Some of these use atmospheric reanalysis data 
provided by the European Centre for Medium-Range 
Weather Forecasts (ECM WF) and National Centers for Envi- 
ronmental Prediction to establish global non-meteorological 
empirical ZTD models such as SHAO-H and IGGTrop??: !°1, 
Some use the Global Geodetic Observing System (GGOS) 
Atmosphere data and International GNSS Service (IGS) 
ZTD data for modeling ZTD directly, such as the GZTD 
and GGZTD models!!!.!7], Additionally, regional ZTD 
models have been developed for China and surrounding 
areasl!3-15], 

The above-mentioned ZTD models are fully empirical, 
established by global atmospheric reanalysis or ZTD data, 
and their global average accuracy is high. However, the accu- 
racy and applicability of ZTD correction in local areas may 
still not be sufficient. To improve the accuracy and applicabil- 
ity of the ZTD model in Asia, this paper establishes a seasonal 
ZTD model, denoted as SSA (Saastamoinen model with sea- 
sonal corrections), which takes into account the annual and 
semi-annual variations based on the widely used Saasta- 
moinen model, using the ZTD grid product with a spatial reso- 
lution of 1° x 1° and a temporal resolution of 6 h released 
by the GGOS Atmosphere. This SSA model is combined 
with a new version of the global pressure-temperature 
model, GPT3!!°l, to form a predictive model for ZTD correla- 
tion called SSA+GPT3, which does not require measured 
meteorological parameters. Using high-precision ZTD data 
from 66 GNSS tracking stations in the Asian region, provided 
by GGOS Atmosphere and IGS as references, the ZTD correc- 
tion accuracy of the SSA model is evaluated for the Asian 
region, and the spatial and temporal distribution characteris- 
tics of the ZTD correction errors are analyzed. 


2. DATA AND METHODS 


2.1. Data Source 


GGOS Atmosphere provides ZTD grid data with a spatial 
resolution of 1° x 1° and a temporal resolution of 6 h based 
on the ECMWF atmospheric reanalysis data from 1980 to 
the present, available from Vienna Mapping Functions 
Open Access Datal!®l, The GGOS Atmosphere products not 
only have a very high spatial resolution, but are also of high 
quality and therefore can be used as a standard dataset for 
ZTD research and applications!!7. 


2.2. Seasonal Effects on the Saastamoinen Model 


The high-precision ZTD provided by the GGOS Atmo- 
sphere from 2015 to 2017 for 66 IGS GNSS tracking stations 
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distributed in the Asian region is used as a reference value 
to calculate the daily averaged bias of the Saastamoinen 
model. Spectral analysis of the bias from 2015 to 2017 is con- 
ducted separately, and the daily averaged bias of these IGS 
stations in each latitude range shows significant annual and 
semi-annual seasonal variations. These variations and their 
spectral analysis results for 2015-2017 at the eastern 
YSSK, southern CUUT, western URUM and northern 
PETS stations in the Asian region are given in Fig. 1. Varia- 
tions at other stations are similar. As can be seen from 
Fig. 1, the daily bias of the Saastamoinen model has annual 
and semi-annual seasonal effects, and these can be approxi- 
mated using the cosine function, 


Bias =BidSmean + A, COS (2x 


a 


doy —d, E 
365.25 


(1) 


Ai cos|2n—>—— 
105 "782.625 


in which the bias and day of year (doy) are known quanti- 
ties; BidSmean is the annual mean value of bias; A; and A2 
are the amplitudes of the annual and semi-annual variations, 
respectively; and dı and dz are the phases of the annual and 
semi-annual variations, respectively. These amplitudes and 
phases can be determined with the nonlinear least squares 
method. The spatial and temporal distribution characteristics 
of the amplitudes and phases are not discussed or analyzed 
here, owing to space limitations. 


2.3. The SSA Model 


Considering the seasonal effects on the Saastamoinen 
model, this paper proposes a new seasonal ZTD model by 
adding the annual and semi-annual terms into the Saasta- 
moinen model, referred to as SSA, 


ZT D =0.002277x 


P+ [00s + 5). 


+ Bias, (2) 
1 — 0.00266 cos(2y) — 0.00028h 


in which P, T, e, and h are the surface pressure in hPa, 
temperature in Kelvin, water vapor pressure in hPa, station lat- 
itude in radians and height in km, respectively. 

To overcome the dependence of the SSA model on mea- 
sured meteorological parameters and improve the applicabil- 
ity of the model, the GPT3 model is introduced to provide rel- 
atively accurate pressure, temperature and water vapor pres- 
sure values. To accomplish this, we construct a refined 
SSA+GPT3 model using the following process. 

Step 1. Calculating the pressure, temperature and water 
vapor pressure parameters of the four horizontal grid points 
nearest to the station, using the GPT3 model, based on the sta- 
tion latitude, longitude, altitude and time information, while 
finding the ZTD, of the four points adjacent to the station in 
the GGOS Atmosphere grid file. 

Step 2. Using equation 3 to make altitude corrections to 
the meteorological parameters of the grid points, we calculate 
the ZT D(h) of each of the four grid points relative to the sta- 
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Fig. 1. Time-series of the daily averaged bias and frequency spectrum analysis between 2015 and 2017 at four stations. (A) YSSK 
(47.03°N, 142.717°E, 91.719 m). (B) CUUT (13.736°N, 100.534°E, 74.699 m). (C) URUM (43.808°N, 87.601°E, 859.352 m). (D) 
PETS (53.023°N, 158.65°E, 102.604 m). (E) YSSK (47.03°N, 142.717°E, 91.719 m). (F) CUUT(13.736°N, 100.534°E, 74.699 m). 


(G) URUM (43.808°N, 87.601°E, 859.352 m). (H) PETS (53.023°N, 


tion altitude according to ZT D(h) = ZT Do - e®", in which 8 
is the height imputation factor, B =—1.3137x10-+. We 
then use bilinear interpolation to calculate the meteorological 
parameters and ZTD at the station, expressed asl!®] 


T=T)+4dT xdh 
P=Po x exp {8m XdMye/[RgXTox(1+0.60770)]xdh} , (3) 
e = eo X (100 x P/Po)**! 


in whichPo, To and eo represent the pressure, temperature 
and water vapor pressure at the grid point surrounding the 
sites, respectively. Q is the specific humidity in kg/kg, 4 is 
the water vapor descent factor, dT is the temperature lapse 
rate in degrees/m, dh is the difference of the site height and 
grid height, gm is the mean acceleration due to gravity, 
&m = 9.80665 m/s”, Rg, is the gas constant, R, = 8.3143 
J/(K-mol), My represents the molar mass of dry air, 
My = 0.028965 kg/mol. 

Step 3. Calculating the ZTD by substituting the meteoro- 
logical parameters at the stations obtained in step 2 into the 
Saastamoinen model and treating them as true values to 
obtain the bias of the Saastamoinen model. 

Step 4. Periodic fitting of the ZTD bias is performed to 
determine the amplitude and phase of the annual and semi- 
annual variations. According to the above model construction 
process, the combined SSA+GPT3 model can perform a 
ZTD calculation at any position. 


3. VALIDATION OF THE MODEL 


3.1. Research Area 


66 IGS GNSS tracking stations in the Asian region are 
selected to analyze and evaluate the effectiveness of the 
model. The range of the research area is 10°N—63°N and 


158.65°E, 102.604 m). 


66°E-159°E. The geographical distribution of the 66 IGS sta- 
tions in the Asian region is shown in Fig. 2. The ZTD products 
from 2015 to 2017 provided by GGOS Atmosphere are 
used as reference values to analyze the spatial and temporal 
distribution characteristics of the SSA model errors. Final 
ZTD data in 2018 from the IGS are used to evaluate the accu- 
racy of the ZTD prediction values of the SSA model, and 
the bias and RMS error are taken as error measurements. 


3.2. Temporal Distribution of Bias and RMS 


Saeed 


To analyze the day-by-day variation of bias and RMS 
of SSA and Saastamoinen models, the daily averaged bias 
and RMS of each IGS station are calculated by day of year. 
The daily averaged bias and RMS variations in 2015-2017 
calculated from SSA+GPT3 and Saastamoinen+GPT3 (abbre- 
viated to SA+GPT3) for the eastern YSSK, southern 
CUUT, western URUM and northern PETS stations in the 
Asian region are shown in Fig. 3. Variations at other stations 
are similar and therefore are not given. 

The daily averaged bias and RMS of the Saastamoinen 
model for the four IGS stations show clear seasonal varia- 
tions, i.e., small in winter and large in summer. The reason 
for these variations is the decrease in wet delay due to dry 
air in winter, which causes a decrease in the error in ZTD esti- 
mates, whereas the wet delay increases because of humid 
air in summer, and hence the error also increases. The bias 
and RMS at the western URM station are relatively small 
compared with those at other stations in Asia, whereas the 
bias and RMS at the eastern YSSK station are larger in sum- 
mer, owing to the effect of the marine climate causing drastic 
variations of water vapor. In contrast to other stations, the 
bias and RMS change rapidly at the southern CUUT station, 
which is located at a low latitude and influenced by the tropi- 
cal climate and marine climate. These results illustrate that 
Saastamoinen model can be improved if these seasonal varia- 


Variations of daily averaged bias and RMS 
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Fig. 2. Distribution of the 66 IGS tracking stations across Asia. 
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Fig. 3. Variations of the daily averaged bias and RMS between 2015 and 2017 at four stations. (A) YSSK (47.03°N,142.717°E, 
91.719 m). (B) CUUT (13.736°N,100.534°E, 74.699 m). (C) URUM (43.808°N,87.601°E, 859.352 m). (D) PETS (53.023°N,158.65°E, 
102.604 m). (E) YSSK (47.03°N, 142.717°E, 91.719 m). (F) CUUT (13.736°N,100.534°E, 74.699 m). (G) URUM (43.808°N,87.601°E, 


859.352 m). (H) PETS (53.023°N,158.65°E, 102.604 m). 


tions are taken into account in ZTD estimation. The bias 
and RMS of the SSA are not only smaller than those of the 
Saastamoinen model, but also more stable, demonstrating 
that the introduction of annual and semi-annual terms into 
ZTD estimates can yield noticeable improvements. 


Dedede 


To analyze the month-to-month variations of the bias 
and RMS of the two models, the monthly averaged bias and 
RMS of 66 IGS stations in the Asian region from 2015 to 
2017 are calculated, as shown in Fig. 4. The bias of the Saas- 
tamoinen model is generally positive, indicating that this 
model usually underestimates the ZTD in the Asian region 
and has a systematic bias. By comparison the monthly aver- 
aged bias and RMS from the SSA model are relatively 
small and stable, with no clear periodic variation in bias, 


Variations of monthly averaged bias and RMS 
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while periodic variation in RMS is substantially reduced. 
The bias and RMS of the SSA do not show significant fluctua- 
tions even in summer and autumn when water vapor levels 
change dramatically, indicating that the SSA model, with sea- 
sonal terms, has a higher accuracy than the conventional Saas- 
tamoinen model, especially in summer and autumn. 


3.3. Spatial Distributions of Bias and RMS 


The 3-year averaged bias and RMS statistics from 2015 
to 2017 are calculated and shown for 66 IGS stations in the 
Asian region in Fig. 5. Here, the bias and RMS of the Saasta- 
moinen model are relatively small in the high-latitude 
region of western and northern Asia, but are larger in the 
southern low-latitude region characterized by oceanic and 
tropical climates. The bias and RMS are also relatively 
large in the eastern coastal region, which is influenced by 


an oceanic climate with drastic water vapor changes. Further- 
more, the bias of the Saastamoinen model is positive at 
most stations in the Asian region, again indicating that this 
model often underestimates the actual value of ZTD. The 
bias and RMS of the SSA model are smaller than those of 
the Saastamoinen model. Specifically, the bias is close to 
zero, indicating that the SSA model, with seasonal terms, 
has little systematic bias. The mean values of the bias and 
RMS of all 66 IGS stations are 1.91 cm and 5.18 cm, respec- 
tively, between 2015 and 2017 for the Saastamoinen model, 
while the mean value of the RMS is 4.2 cm for the SSA, 
with an overall reduction of 19% relative to the Saasta- 
moinen. 


3.3.1. Variations of bias and RMS with height 


To analyze the characteristics of height distribution of 
the bias and RMS of the two models, the altitude of 66 sta- 
tionswas dividedintobrackets of<500m,>500—1000m, 1500- 
2 000 m and >2 000 m. The bias and RMS at each height 
range from 2015 to 2017 are shown in Fig. 6, in which the 
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bias from the SSA is small enough to be insignificant, and 
has been omitted. The RMS of both models shows a notice- 
able decreasing trend with increasing station height. The 
RMS of the SSA is smaller than that of the Saastamoinen 
model in each height range. Also noteworthy is the fact that 
the bias of the Saastamoinen model is positive from 0 to 
2000 m and becomes larger with increasing height. How- 
ever, the bias shows a sudden increase when the altitude 
exceeds 2000 m. One possible explanation is that there are 
few IGS stations in this height range. Furthermore, the accu- 
racy of both models is better at high altitudes. This is 
because the integration interval of tropospheric refractivity 
is shortened at high altitudes, with drier air causing reduced 
wet delay at high-altitude stations, such as URUM and 
LHAZ in China. In all latitude ranges, the accuracy of the 
SSA is improved to different degrees compared with the Saas- 
tamoinen model. Overall, the improvement is most significant 
in low-altitude regions. Compared with the Saastamoinen 
model, the RMS of the SSA model is reduced by 20.31%, 
11.25%, 4.01%, and 11.05% in the height ranges of 
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Fig. 4. Monthly averaged bias and RMS from the SSA and Saastamoinen models. (A) Monthly averaged bias during 2015-2017. 


(B) Monthly averaged RMS during 2015-2017. 
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Fig. 5. Distribution of the bias and RMS between 2015 and 2017 from the SSA and Saastamoinen models. (A) Bias distribution of 
Saastamoinen. (B) RMS distribution of Saastamoinen. (C) Bias distribution of SSA. (D) RMS distribution of SSA. 
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Fig. 6. Variations of the bias and RMS between 2015 and 2017 from the SSA and Saastamoinen models with height. (A) Bias 


during 2015-2017. (B) RMS during 2015-2017. 


<500 m, >500-1000 m, 1500-2 000 m, and >2 000 m, 
respectively. 


B32: 


To analyze the latitude distribution characteristics of 
the bias and RMS of both models, the latitudes of 66 stations 
were divided into 16°N-30°N, 30°N-45°N, and 45°N-— 
63°N regions. The averaged bias and RMS of 66 stations in 
each latitude range from 2015 to 2017 are shown in Fig. 7, 
in which the bias of the SSA is very low and has been omit- 
ted. The RMS shows a decreasing trend with increasing lati- 
tude, and the accuracy of the SSA is superior to that of the 
Saastamoinen model in all latitude ranges. The bias of the 
Saastamoinen model also shows a decreasing trend with 
increasing station latitude. Additionally, the accuracy of the 
SSA is substantially better at the high-latitude area above 
30°N compared with that of the Saastamoinen. This is 
because the ZTD estimation error is larger in the lower lati- 
tudes of Asia because of the oceanic and tropical climates, 
where water vapor changes are more complex. The accuracy 
of the SSA model improves to different degrees in all latitude 
ranges compared with that of the Saastamoinen model, and 


Variations of bias and RMS with latitude 


A 3.5 


3.30 mm 16°N-30°N 
EE 30°N-45°N 
3.0+ E 45°N-60°N 
2.54 
£2.01 
2 
3 
515l 
1.15 
Lol 100 
j l B 
0.0 
SA+GPT3 
Model 


this enhancement is more obvious in regions below 30°N. 
In contrast with the Saastamoinen model, the RMS of the 
SSA model decreases by 30%, 9.55% and 10.3% in the lati- 
tude ranges from 16°N—30°N, 30°N—45°N and 45°N—63°N, 
respectively. 


3.4. Prediction Assessment of the SSA 


To further assess the accuracy of the SSA model, the 
model is tested with the ZTD grid products released by 
GGOS Atmosphere between 2015 and 2017. This SSA 
model is then employed to predict the ZTD of 66 IGS stations 
in the Asian region in 2018. The final ZTD data from the 
IGS are then taken as the real values for model assessment. 
Fig. 8 shows the spatial distribution of the yearly averaged 
bias and RMS of the ZTD predictions, and the bias and 
RMS statistical measures are given in Table 1. 

The yearly averaged bias and RMS from the Saasta- 
moinen model are smaller in the western region and high-lati- 
tude region of Asia, and larger in the eastern and low-latitude 
regions, with a maximum bias of 6.07 cm. Compared with 
the Saastamoinen model, the yearly averaged bias of the 
SSA is smaller in the whole Asian region with a maximum 
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Fig. 7. Variations of the bias and RMS between 2015 and 2017 from the SSA and Saastamoinen models with latitude. (A) Bias 


during 2015-2017. (B) RMS during 2015-2017. 
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Fig. 8. Distribution of the yearly bias and RMS of the ZTD predictions from the SSA and Saastamoinen models. (A) Bias 
distribution of Saastamoinen. (B) RMS distribution of Saastamoinen. (C) Bias distribution of SSA. (D) RMS distribution of SSA . 


Table 1. Statistics of the yearly bias and RMS of the ZTD 
predictions over Asia 


na Saastamoinen SSA 
Statistics 
Mean Range Mean Range 
Bias/cm 1.45 [-2.71, 6.07] —0.38  [=2.43, 1.19] 
RMS/em 5.16 [1.74, 8.44] 4.43 [1.76, 6.83] 


value of 2.43 cm. The RMS is approximately 5 cm in north- 
western Asia and larger in other regions for the Saastamoinen 
model, whereas the value for the SSA is under 5 cm in 
regions with latitudes above 40°N and larger in lower lati- 
tudes. Additionally, the enhancement of the SSA is most obvi- 
ous in lower-latitude regions. The RMS of both models 
becomes smaller with increasing latitude. The SSA model 
reduces the bias and RMS by 73.79% and 14.15%, respec- 
tively, compared with the Saastamoinen model. 

To further analyze seasonal variations of the bias and 
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RMS, the monthly averaged bias and RMS of 66 IGS stations 
in 2018 is shown in Fig. 9. The monthly averaged bias and 
RMS of the Saastamoinen model show obvious seasonal vari- 
ations, large in summer and small in winter, whereas the 
RMS of the SSA model changes slowly in each month and 
does not show large fluctuations even in the summer and 
autumn. The bias of the SSA, shown in Fig. 9A, is signifi- 
cantly smaller than that of the Saastamoinen model except 
in January, March, and October. This is because the bias of 
the SSA model has more negative values in these months 
than the Saastamoinen model, leading to a large calculated 
monthly averaged bias. However, although the monthly aver- 
aged RMS of both models is comparable in spring and win- 
ter, the monthly averaged RMS of the SSA in autumn and 
summer is significantly lower than that of the Saastamoinen 
model, and the RMS is reduced by 36.36% at maximum. 
Fig. 8, Fig. 9 and Table | not only show that the ZTD predic- 
tion accuracy of the SSA in the Asian region is better than 
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Fig. 9. Seasonal variations in the bias and RMS of the ZTD predictions from the SA+GPT3 and SSA+GPT3 models over Asia. 
CA) Monthly averaged bias in 2018. (B) Monthly averaged RMS in 2018. 
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that of the Saastamoinen model, but also verify that the spatial 
and temporal variation characteristics of the prediction 
errors are generalized. 


4. CONCLUSIONS 


The widely used Saastamoinen model requires in-situ 
meteorological parameters, and it is also affected by seasonal 
changes in environmental conditions. In this study, the bias 
of the Saastamoinen model shows clear seasonal variations, 
and the spectrum analysis indicates that these variations are 
mainly dominated by annual and semi-annual variations. To 
improve the Saastamoinen model we propose an enhanced 
model, dubbed the SSA model, to compensate for seasonal 
annual and semi-annual variations, using meteorological 
parameters from the empirical GPT3 model. The results 
demonstrate that the SSA, with annual and semi-annual 
terms, has a high ZTD prediction accuracy in the Asian 
region with a yearly averaged RMS of approximately 4.5 
cm. Both the bias and RMS of the SSA model are smaller 
and more stable than those of the Saastamoinen model, 
even in summer and autumn. The enhancement in summer 
can reach 36.36% in contrast to the Saastamoinen model. 
The extent of this improvement is also related to the geo- 
graphic location; the improvement is very noticeable in the 
low-latitude and low-altitude areas of Asia. Moreover, the 
accuracy enhances with an increase in latitude or altitude. 
Because the SSA model is constructed based on 1° x 1° 
high-resolution GGOS Atmosphere data products, and the 
required meteorological parameters can be provided by 
GPT3, we conclude that the combined SSA+GPT3 model is 
a potentially very effective method for tropospheric delay cor- 
rections in high-precision GNSS applications over Asia, 
and can be used to predict tropospheric delays in real-time 
high-precision GNSS positioning, navigation and timing. 
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